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Hydropower  is  not  only  a  renewable  and  sustainable  energy  source,  but  its  flexibility  and  storage 
capacity  also  makes  it  possible  to  improve  grid  stability  and  to  support  the  deployment  of  other 
intermittent  renewable  energy  sources  such  as  wind  and  solar  power.  As  a  result,  a  renewed  interest  in 
pumped-hydro  energy  storage  plants  (PHES)  and  a  huge  demand  for  the  rehabilitation  of  old  small 
hydropower  plants  are  emerging  globally. 

As  regards  PHES,  advances  in  turbine  design  are  required  to  increase  plant  performance  and  flexibility 
and  new  strategies  for  optimizing  storage  capacity  and  for  maximizing  plant  profitability  in  the 
deregulated  energy  market  have  to  be  developed. 

During  the  upgrading  of  old  small  hydropower  plants,  the  main  challenges  to  be  faced  are  the  design 
of  new  runners,  that  had  to  match  the  existing  stationary  parts,  and  the  development  of  optimal  sizing 
and  management  strategies  to  increase  their  economic  appeal. 

This  paper  traces  an  overview  of  the  prospects  of  pumped-hydro  energy  storage  plants  and  small 
hydro  power  plants  in  the  light  of  sustainable  development.  Advances  and  future  challenges  in  both 
turbine  design  and  plant  planning  and  management  are  proposed.  PHES  and  hybrid  wind/solar-PHES  are 
illustrated  and  discussed,  as  well  as  the  limits  and  peculiarities  of  the  new  design  strategies,  based  on 
computational  fluid  dynamics,  for  both  PHES  and  small  hydropower  plants. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 
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In  the  last  decades,  world  electrical  energy  consumption  has 
significantly  increased  with  a  share  that  has  reached  17.7%  in  2010 
and  is  predicted  to  double  by  2025  [1]. 

The  increasing  concern  about  environmental  aspects  has 
favored  a  corresponding  rapid  growth  of  the  deployment  of  the 
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renewable  energy  sources  aimed  at  the  progressive  reduction  of 
fossil  fuel  exploitation  and  dependence.  In  such  a  context  hydro- 
power  is  undoubtedly  one  of  the  most  mature  technologies  with 
an  electricity  production  of  about  3500  TWh  in  2010  (16.3%  of  the 
world's  electricity),  greater  than  that  of  the  other  renewable 
sources  combined  (3.6%),  but  much  smaller  than  that  of  the  fossil 
fuel  plants  (67.2%)  [2]  (Fig.  1). 


In  the  EU-27  the  hydropower  generation  was  323  TWh  in  2010 
(9.8%  of  the  European  electricity)  and  around  60%  of  electricity 
generation  from  renewable  sources  [3], 

Even  though  today  hydropower  plays  a  key  role  in  the  green 
energy  production,  avoiding  the  combustion  of  4.4  million  barrels 
of  oil  equivalent  daily,  only  33%  of  potential  hydro  resources  has 
been  developed  and  the  remaining  technical  potential  is  estimated 
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Fig.  1.  World  electricity  generation  by  fuel. 
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Fig.  2.  World  hydropower  technical  potential  [2], 
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to  be  very  high  (14,576  TWh/year)  [2]  (Fig.  2).  The  highest 
percentage  of  undeveloped  potential  is  located  in  Africa  (92%), 
followed  by  Asia  (80%),  Oceania  (80%)  and  Latin  America  (74%), 
even  though  this  region  is  also  characterized  by  two  of  the  top  ten 
hydropower  producers  in  2010  (Table  1)  and  by  several  countries 
with  a  very  large  share  of  electricity  generation  obtained  by 
hydropower  (Table  2)  [2], 

Besides  the  positive  effects  on  climate  mitigation,  hydropower 
also  presents  other  considerable  advantages:  it  promotes  price 
stability  because,  unlike  fuel  and  natural  gas,  it  is  not  subject  to 
market  fluctuations;  it  reduces  environment  vulnerability  to 
floods;  it  contributes  to  fresh  water  storage  for  drinking  and 
irrigation  exploitations;  it  makes  a  significant  contribution  to 
development  by  bringing  electricity,  roads,  industry  and  com¬ 
merce  to  communities  which  can  benefit  future  generations  as 
hydropower  projects  are  long-term  investments  with  an  average 
life  span  of  50-100  years. 

Despites  numerous  examples  of  excellent,  sustainable  and  safe 
exploitation  of  the  water  resources  by  hydropower,  large  hydro 
projects,  which  include  a  dam  and  a  reservoir,  encountered  a 
substantial  opposition  in  the  latter  end  of  the  last  century  because 
of  their  environmental  and  social  implications  (landscape,  wildlife, 
biodiversity,  population  settlement,  health  and  water  quality,  etc.). 
This  opposition  was  one  of  the  main  factors  of  the  slowdown  in 
the  hydroelectricity  generation  between  the  late  1990s  and  the 
early  2000s  [2,3], 

However,  eliminating  large  hydropower  projects  from  renew¬ 
able  energy  programs  will  not  reduce  the  power  demand,  which 
will  be  partially  satisfied  by  thermal  plants,  thereby  increasing  the 
global  level  of  greenhouse  gas  emissions.  For  instance,  it  was 
demonstrated  that  two  pumped-hydro  energy  storage  units  com¬ 
bined  with  a  thermal  generation  unit  make  it  possible  to  reduce 
the  excess  emissions  of  the  thermal  unit  by  60%  [4], 

In  such  a  context,  a  renewed  interest  in  large  pumped-hydro 
energy  storage  plants  and  a  huge  demand  for  the  rehabilitation 
and  repowering  of  old  small  and  micro-hydro  plants  are  globally 
emerging  both  due  to  further  increases  in  the  corresponding  share 
of  renewable  electricity  production  and  due  to  the  support,  in 
terms  of  storage  capacity,  of  a  wider  exploitation  of  other  renew¬ 
able  energy  sources  such  as  wind  and  solar  power. 


Table  1 

Top  ten  hydropower  producers  in  2010  [2], 


Country 

Hydro  electricy  [TWh] 

Share  of  electricity  generation  [%] 

China 

694 

14.8 

Brazil 

403 

80.2 

Canada 

376 

62.0 

United  States 

328 

7.6 

Russia 

165 

15.7 

India 

132 

13.1 

Norway 

122 

95.3 

Japan 

85 

7.8 

Venezuela 

84 

68.0 

Sweden 

67 

42.2 

Even  though  these  renewable  energy  sources  have  largely 
increased  in  the  last  years,  reaching  in  2010  a  share  of  3.6%  of 
the  world  electricity  generation  [1],  this  share  would  need  to  be 
greatly  increased  in  order  to  significantly  contribute  to  the  reduc¬ 
tion  of  greenhouse  gas  emissions  (GHG).  The  European  Union 
policy  has  made  great  efforts  towards  developing  these  resources 
in  response  to  environmental  concerns,  but  their  unpredictable 
and  intermittent  characteristics  have  restrained  their  deployment 
because  of  the  negative  impact  on  power  system  security,  stability, 
reliability  and  efficiency  [5-8]. 

At  present,  the  electricity  grid  is  highly  centralized  with  a 
complex  system  of  energy  production-transmission  characterized 
by  a  long  distance  between  power  plants  and  end-users  and  by  a 
limited  use  of  storage,  whose  installed  capacity  was  about 
127.9  GW  in  2010  (2.5%  of  the  world  installed  capacity)  [9].  To 
ensure  the  security  of  the  power  system,  a  continuous  balance 
between  demand  and  supply  should  be  guaranteed  and  this 
actually  limits  penetration  in  the  grid  of  intermittent  renewable 
energy  sources,  whose  energy  production  is  fluctuating,  unpre¬ 
dictable  and  delocalized.  For  example,  as  regards  the  potential 
wind  energy  penetration  in  the  electricity  grid,  an  instantaneous 
increase  up  to  20%  of  the  total  energy  production  was  estimated  to 
be  feasible  without  technical  hitches  for  the  grid  stability  [10], 
However,  further  increases  in  wind  energy  generation  (up  to  80%) 
are  feasible  only  in  isolated  grids  (smaller  than  10  MW),  whereas 
greater  electrical  grids  require  energy  storage  to  accept  a  wind 
energy  penetration  greater  than  20%  [11], 

The  development  of  a  significant  energy  storage  capacity  is, 
therefore,  a  necessary  solution  to  favour  the  deployment  of  the 
renewable  energy  sources  not  only  in  isolated  grids,  but  also  in 
interconnected  grid  systems,  as  demonstrated  by  several  analyses 
carried  out  on  a  national  scale  [12-17]. 

For  this  reason,  the  European  Union  is  carrying  out  a  Climate 
and  Energy  policy,  defined  in  the  Strategic  European  Technology 
Plan  (SET-Plan)  [3],  one  goal  of  which  is  to  study  more  in  depth  the 
benefits  of  storage  applications.  In  such  a  context,  several  studies 
[9,18-23]  have  been  carried  out  to  analyse  the  current  status  of  the 
wide  range  of  available  technologies  (mechanical,  electromagnetic, 
chemical,  thermal)  in  terms  of  technology  maturity,  efficiency, 
energy  storage  capacity,  power  discharged  capacity,  application 
size,  cost  of  investment,  life  time  and  environmental  impact. 

All  these  analyses  identify  Pumped  Hydro  Energy  Storage  (PHES) 
as  being  the  most  cost-efficient  large-scale  storage  technology 
currently  available,  with  an  efficiency  range  of  75-85%  and  compe¬ 
titive  costs  (1500-4300  $/l<W,  250-430  $/kWh)  [23],  In  Europe,  this 
technology  represents  99%  of  the  on-grid  electricity  storage  [9]  with 
more  than  7400  MW  of  new  PHES  plants  proposed  and  a  total 
investment  cost  of  over  €  6  billion  [24],  In  spite  of  this  boost  provided 
by  the  increasing  need  for  storage  capacity,  one  of  the  major  limits  for 
the  further  development  of  PHES  is  the  lack  of  suitable  locations  for 
the  construction  of  new  facilities  [19,25-27].  To  overcome  this 
problem,  a  program,  able  to  identify  areas  that  could  be  quite  easily 
modified  in  order  to  construct  the  reservoirs  of  the  PHES,  was 
developed  [28]  and  the  possibility  of  excavating  underground  reser¬ 
voirs  was  considered  [29,30],  However,  in  the  long-term  perspective 


Table  2 

Countries  with  more  than  a  half  of  their  electricity  generation  from  hydropower  in  2010  [2], 


Share  of  hydropower 

Countries 

Hydropower  Generation 

-100% 

Albania,  DR  of  Congo,  Mozambique,  Nepal,  Paraguay,  Tajikistan,  Zambia 

54  TWh 

>90% 

Norway 

126  TWh 

>80% 

Brazil,  Ethiopia,  Georgia,  Kyrgyzstan,  Namibia 

403  TWh 

>70% 

Angola,  Columbia,  Costa  Rica,  Ghana,  Myanmar,  Venezuela 

77  TWh 

>60% 

Austria,  Cameroon,  Canada,  Congo,  Iceland,  Latvia,  Peru,  Tanzania,  Togo 

351  TWh 

>50% 

Croatia,  Ecuador,  Gabon,  DPR  of  Korea,  New  Zealand,  Switzerland,  Uruguay,  Zimbabwe 

36  TWh 
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to  fulfil  the  increasing  storage  capacity  need,  besides  the  installation 
of  new  PHES  sites,  it  would  be  necessary  to  adapt  and  exploit  the 
existing  hydropower  plants,  as  was  achieved  in  France  during  the 
1970s  and  1980s  to  support  the  nuclear  power  reactors  [22]  or  as 
recently  proposed  in  Greece  to  support  RES  penetration  [12],  To  reach 
this  aim,  the  existing  hydropower  and  PHES  plants  should  not  only 
modify  their  operation  strategies  in  order  to  maximize  the  revenue 
from  the  day-ahead  market  and  from  the  regulation  services,  but  also 
optimize  turbine  performance  by  means  of  innovative  design  criteria 
so  as  to  increase  the  corresponding  storage  efficiency  and  to  achieve 
the  required  greater  flexibility.  The  resulting  storage  capacity  will  also 
favour  a  different  economical  approach  to  investment  in  RES  produc¬ 
tion  based  not  only  on  the  incentive  mechanisms,  but  also  on  the 
techno-economic  optimization  of  the  plant  operation.  This  would  be 
aimed  at  maximizing  returns  by  storing  the  amount  of  green 
production  during  the  low  demand  periods  and  selling  it  to  the 
system  during  the  peak  demand  periods,  thereby  providing  the 
required  grid  stability  by  means  of  a  fast  response  time. 

As  regards  small  hydropower  plants  (SHP),  the  main  barriers  to 
the  development  of  new  plants  are  represented  by  the  environ¬ 
mental  requirements,  in  particular  those  defined  by  the  Environ¬ 
mental  Impact  Assessments  (EIA)  and  by  the  Water  Framework 
Directive  (WFD).  These  requirements,  generally  extremely  restric¬ 
tive,  do  not  properly  consider  the  benefits  of  the  SHP  renewable 
energy  production  and  lead  to  a  limitation  of  the  hours  of 
production  and  to  an  increase  of  the  investment  costs.  So,  to 
safeguard  the  profitability  of  the  SHP  plants  and  to  obtain  a  further 
increase  of  the  green  energy  production,  the  most  cost-effective 
solutions  are  represented  on  the  technical  side  by  an  increase  of 
the  plant  efficiency  by  runner  upgrade,  and  on  the  economical  side 
by  an  optimal  sizing  of  the  plant  accompanied  by  the  definition  of 
optimal  management  strategies. 

The  paper  is  structured  as  follows:  Section  2  focuses  on 
Pumped  Hydro  Energy  Storage  Plants.  Following  a  brief  introduc¬ 
tion  to  technology  advances  and  challenges,  the  latest  trends  in 
the  definition  of  optimal  operating  strategies  are  reviewed: 
standard  configurations  as  well  as  innovative  hybrid  wind-solar/ 
pumped  hydro  configurations  are  considered.  Section  3  deals  with 
small  hydro  power  plants  and,  after  an  overview  of  the  main 
characteristics  of  hydraulic  turbines,  reviews  advances  and  trends 
in  optimal  sizing  and  management  of  these  type  of  plants.  Section 
4  focuses  on  the  development  of  the  computational  fluid- 
dynamics  for  turbomachinery  design  and  in  particular  on  its  key 
role  as  an  innovative  and  effective  design  tool  in  the  rehabilitation 
of  small  hydroelectric  plants  and  in  the  development  of  the  new 
generation  of  pumped-hydro  energy  storage  plants. 


2.  Pumped  hydro  energy  storage  plants  (PHES) 

A  hydro  pumped  energy  storage  plant  converts  grid- 
interconnected  electricity  to  hydraulic  potential  energy  (so-called 
“charging”),  by  pumping  the  water  from  a  lower  reservoir  to  an 
upper  one  during  the  off-peak  periods,  and  then  converting  it  back 
during  the  peak  periods  (“discharging”)  by  exploiting  the  available 
hydraulic  potential  energy  between  the  reservoirs  like  a  conven¬ 
tional  hydropower  plant. 

Pumped  energy  storage  plants  are  generally  subdivided  in 
“closed-loop”  (or  "off-stream”)  plants,  when  the  discharging  phase 
relies  entirely  on  the  water  previously  pumped  to  the  upper 
reservoir,  and  in  “pump-back”  ones,  when  the  discharging  phase 
relies  on  a  combination  of  pumped  water  and  natural  inflow 
(Fig.  3). 

These  plants  require  very  specific  site  conditions  to  be  feasible 
and  viable,  among  which  proper  ground  conformation,  differe¬ 
nce  in  elevation  between  the  reservoirs  and  water  availability. 


a 


b 


Fig.  3.  Closed-loop  (a)  and  pump-back  (b)  pumped  energy  storage  plants. 


For  these  reasons,  the  earliest  PHES  plants  were  built  in  the  Alpine 
regions  of  Switzerland  and  Austria  whose  ground  conformation 
together  with  the  presence  of  hydro  resources  was  suitable 
for  PHES. 

The  first  PHES  plants,  owned  by  state  utilities,  were  built  to 
supply  energy  during  the  peak  periods,  allowing  the  base-load 
power  plants  to  operate  at  high  efficiency,  and  to  provide  balan¬ 
cing,  frequency  stability  and  black  starts.  The  period  from  the 
1960s  to  the  late  1980s  was  characterized  by  a  significant 
development  of  these  plants  mainly  due  to  the  corresponding 
deployment  of  nuclear  power  plants  whose  great  inertia  was 
compensated  by  the  PHES  flexibility  [24],  In  the  1990s  the  reduced 
growth  of  the  nuclear  plants  together  with  the  increasing  difficult 
identification  of  suitable  locations  significantly  limited  the  further 
development  of  new  PHES  plants.  For  this  reason,  in  2006  the 
average  percentage  of  PHES  installed  capacity  was  only  about  6% 
of  the  full  generation  installed  capacity  in  the  majority  of  the 
world  countries  with  the  exception  of  Luxemburg  (67%)  [24],  This 
percentage  was  greater  than  10%  only  for  those  countries  char¬ 
acterized  by  a  significant  availability  of  hydro  sources  (Croatia)  or 
by  a  significant  percentage  of  installed  nuclear  power  capacity 
(Latvia,  Japan  and  the  Slovak  Republic).  USA  and  Japan  still 
maintained  the  world  highest  installed  PHES  capacity  with 
20815  MW  and  24575  MW  respectively,  whereas  in  the  European 
context  the  largest  number  (23)  of  PHES  plants  are  concentrated  in 
Germany  [24], 

2.3.  Variable-speed  technology 

In  recent  years,  after  the  liberalization  of  the  market,  the 
increasing  interest  in  renewable  energy  sources  has  again  turned 
public  attention  towards  the  PHES  as  a  mature  and  large  scale 
energy  storage  technology  to  support  green  energy  production 
and  to  provide  grid  stability.  In  such  a  context,  several  new 
PHES  have  been  planned  in  Europe  for  a  total  power  capacity 
of  7426  MW  [31]  and  some  of  them  will  adopt  the  variable- 
speed  reversible  pump-turbines  breakthrough  technology,  whose 
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peculiarity  is  to  improve  the  pump-turbine  efficiency  over  a  wider 
range  of  operating  conditions  and  to  improve  the  capability  in  grid 
regulation  of  the  PHES  [32,33]. 

Usually,  pump-turbines  are  optimized  for  a  specific  speed, 
discharge  flow  rate  and  head.  At  this  fixed  speed,  only  limited 
deviations  of  head  and  discharge  are  allowed  since,  at  off-design 
operating  conditions,  the  pump-turbine  efficiency  significantly 
drops.  According  to  the  similarity  theory,  at  different  rotation 
velocities,  the  operating  conditions  characterized  by  the  best 
efficiency  are  different.  The  basic  idea  underlying  variable-speed 
technology  is  to  maintain  high  efficiency  values  in  the  whole 
working  range  of  the  turbine  by  modifying  the  rotation  velocity  of 
the  machine  in  order  to  remain  near  the  best  efficiency  line, 
thereby  maintaining  optimal  or  near-to  optimal  efficiency  values 
in  a  wide  range  of  head  and  flow  rates. 

Fig.  4  reports  a  comparison  of  operating  ranges  in  terms  of  head 
and  discharge  variation  between  a  fixed  and  a  variable-speed 
pump-turbine  with  reference  to  the  pumping  mode.  It  is  clear  that 
the  possibility  of  varying  the  speed  also  significantly  increases  the 
operating  range  of  the  machine,  which  was  limited,  on  the  one 
hand,  by  cavitation  problems  (broken  line  on  the  right)  and  on  the 
other  by  an  instability  operating  area  (broken  line  on  the  left),  a 
typical  behaviour  of  the  pump-turbine  at  part  load  that  should  be 
avoided  during  the  plant  operation  since  it  may  lead  to  self-excited 
vibrations  of  the  hydraulic  system. 

From  an  electrical  point  of  view,  the  possibility  of  speed  variation 
is  obtained  by  means  of  a  power  converter  employing  power 
electronic  to  decouple  motor/generator  from  the  grid  in  terms  of 
reactive  power,  voltage  and  frequency,  that  are  properly  set  indepen¬ 
dently  on  the  different  sides  (power-grid  and  motor/generator). 


The  corresponding  improvements  in  the  design  criteria 
demanded  by  the  variable-speed  technology  do  not  significantly 
increase  the  cost  investment  of  the  power  unit  compared  with  a 
conventional  PHES  plant  [34], 

Two  of  the  largest  PHES  under  construction  in  Europe  that  are 
equipped  with  this  technology  are  Nant  de  Drance  and  Linthal, 
located  in  the  south  west  and  the  north  east  of  Switzerland 
respectively.  In  both  these  power  plants  the  use  of  the  variable 
speed  technology  was  justified  by  the  wide  head  variation:  in  the 
power  station  of  Nant  de  Drance,  equipped  with  six  variable  speed 
units  with  a  unit  output  of  157  MW  (rated  speed=428.6  rpm; 
speed  range=+7%),  the  gross  head  varies  between  250  and 
390  m;  in  the  power  station  of  Linthal,  equipped  with  four  variable 
speed  units  with  a  unit  output  of  250  MW  (rated  speed= 500  rpm; 
speed  range=+6%),  the  gross  head  varies  between  560  and 
724  m  (Table  3).  Table  4  reports  the  project  features  of  another 
important  installation,  that  is  the  PHES  of  Goldisthal,  in  Germany, 
equipped  with  two  fixed  units  and  two  variable-speed  units. 

It  is  interesting  to  highlight  that  in  these  plants  the  design 
priority  was  given  to  flexibility  in  pumping  mode,  maximizing 
the  pump  operating  range,  whereas  in  other  cases,  such  as  the 
PHES  of  Tehri  under-construction  in  India  (four  variable-speed 
units  of  255  MW;  rated  speed  =  230.77  rpm;  speed  range  =  +  6%) 
the  hydraulic  design  was  optimized  to  increase  the  global  plant 
efficiency  (Table  3).  These  different  design  strategies,  aimed  at 
maximizing  the  plant  revenue,  are  due  to  the  different  electricity 
market  regulations  and  in  particular  to  the  existence  (or  not) 
of  a  remunerative  regulation  market,  that  could  significantly 
affected  the  definition  of  the  plant  management  strategies 
(Section  2.2). 


b 


Fig.  4.  Comparison  of  operating  domains  between  fixed  (a)  and  variable-speed  (b)  pump-turbines  (pumping  mode). 
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Table  3 

Project  features  of  the  PHES  of  Nant-de-Drance  and  Linthal. 


Project  features 

Nant  de  Drance  2019 

Linthal  2015 

Tehri  2016 

Rotational  synchronous  speed 

428.6  [rpm] 

500.0  [rpm] 

230.77  [rpm] 

Speed  range 

±7% 

±6% 

±6% 

Head  variation 

250/390  [m] 

560/724  [m] 

830/740  [m] 

Nominal  output  per  unit 

Maximum  pump  discharge  per  unit 

Maximum  turbine  discharge  per  unit 

157  [MW] 

56  [m3/s] 

60  [m3/s] 

250  [MW] 

255  [MW] 

Generator  mode 

175  [MVA] 

280  [MVA] 

278  [MVA] 

Motor  mode 

172  [MW] 

250  [MW] 

Runner  diameter 

6.009  [m] 

4.230  [m] 

Table  4 

Project  features  of  the  PHES  of  Goldisthal. 


Project  features 

Goldisthal 

Rotational  synchronous  speed 

333  [rpm] 

Speed  range  (2  units) 

+4%/- 10% 

Head  variation 

279.2/334.0  [m] 

Nominal  pump  discharge  per  unit 

80  [m3/s] 

Nominal  turbine  discharge  per  unit 

103  [m3/s] 

Nominal  output  per  unit 

265  [MW] 

Runner  diameter 

4.593  [m] 

Even  if  these  plants  represent  the  starting  point  of  a  new 
generation  of  PHES,  the  possibility  in  terms  of  speed  variation  is 
limited  to  about  +  10%  due  to  the  lack  of  full  variable-speed  design 
criteria,  allowing  the  pump-turbine  to  operate  at  a  high  efficiency 
in  a  wide  range  of  rotation  velocities,  and  to  the  lack  of  a  full 
power  converter  technology,  enabling  a  full-range  variable  speed 
turbining  and  pumping  (Tables  3  and  4). 

As  regards  the  power  converter  technology,  significant 
advances  have  been  made  in  recent  years,  as  for  example  the 
100  MVA  variable  speed  frequency  converter  provided  by  ABB  for 
the  PHES  of  Grimsel  2  (Switzerland),  enabling  a  wider  variable 
speed  pumping  (  + 12%)  but  not  turbining.  However,  to  increase 
the  voltage  range  and  to  reduce  cost,  size  and  losses  still  remain 
the  main  challenges  to  face  in  order  to  develop  a  + 100%  power 
converter. 

As  regards  the  mechanical  equipment,  Computational  Fluid 
Dynamics  (CFD)  has  allowed  significant  advances  to  the  under¬ 
standing  of  the  reasons  of  the  unstable  behaviour  of  the  pump- 
turbines  [35-37]  that  prevent  PHES  from  operating  at  low  load 
and  seems  to  be  associated  to  the  development  of  fluid-dynamical 
phenomena  due  to  the  rotor-stator  interaction  (see  Section  4.1). 
However,  innovative  design  criteria  allowing  electricity  production 
in  the  whole  operating  range  (0-100%  of  the  peak  power)  still  have 
to  be  developed  and  certainly  represents  the  future  challenge  for 
the  development  of  a  new  generation  of  pump-turbine. 

2.2.  Optimal  operating  strategies  for  pumped  hydro  energy  storage 
plants 

In  a  liberalized  electricity  market,  pumped  hydro  energy 
storage  represents  a  merchant  unit  whose  revenue  derives  not 
only  from  the  electricity  trading  in  the  day-ahead  market,  but  also 
from  regulation  services  selling  in  the  ancillary  reserve  markets. 

As  regards  electricity  trading,  the  profit  of  the  PHES  is  strictly 
connected  with  selling  hydropower-generated  electricity  at  a  high 
market  clearing  price  (MPC)  in  the  generation  mode  and  with 
purchasing  it  at  a  low  MPC  in  the  pumping  mode.  However,  since 
the  conversion  process  (electricity-mechanical  energy-electricity) 
is  unavoidably  characterized  by  energy  losses,  to  make  revenue, 
the  purchasing  price  should  be  lower  than  the  selling  price  of  a 


value  depending  on  the  plant  conversion  efficiency  (the  so-called 
“round-trip  efficiency”),  that  was  assumed  to  be  at  around  67%  by 
Lu  et  al.  [38]  and  around  75%  by  Castronuovo  et  al.  [39], 

As  for  the  ancillary  markets,  the  participation  of  the  PHES  is 
defined  on  an  hourly  basis  and  is  remunerated  by  an  hourly 
ancillary  service  price  to  which  an  hourly  spot  price  is  added  if 
the  plant  is  called  to  generate  electricity  [38,40], 

The  resulting  profit  of  a  PHES  plant  operating  both  in  the  day- 
ahead  market  and  in  the  regulation  market  is  therefore  repre¬ 
sented  by 

k  —  k generation  "P  ^ regulation  ~  f- pumping  ~  (-0&M  ( 1 ) 

where  R  is  the  plant  profit,  Rgen eration  is  the  energy  revenue 
received  from  the  day-ahead  market,  Regulation  is  the  regulation 
revenue  received  from  the  ancillary  markets,  Cpumping  are  the 
pumping  costs  and  C0&m  are  the  operation  and  maintenance  costs. 
In  a  deregulated  market  with  prices  varying  on  a  hourly  basis,  the 
profit  R[0,t]  of  a  PHES  plant  operating  over  a  period  of  T  hours 
depends  on  the  hourly  prices  Pgen  (€/MWh)  received  for  generat¬ 
ing  Egen  (MWh),  on  the  hourly  spinning  reserve  prices  Prs  (€/MWh) 
received  for  generating  Ers  (MWh)  in  the  regulation  market,  on  the 
hourly  non-synchronous  reserve  prices  Pm  (€/MWh)  received  for 
staying  off-line  without  producing  £m  (MWh),  on  the  hourly 
purchasing  prices  Ppump  (€/MWh)  paid  for  pumping  and  on  the 
operation  and  maintenance  costs,  assumed  to  be  constant 


h 

—  2  EpUmpPpUmp(t)  —  Co&m  (2) 

i=  l 

where  tg,  trs,  tm  and  tp  are  the  number  of  hours  spent  by  the  plant 
in  generating  mode,  in  spinning  regulation  mode,  in  non- 
synchronous  regulation  mode  and  in  pumping  mode  respectively. 

Since  the  price  volatility  of  the  electricity  market  does  not 
necessarily  follow  the  load  trend,  proper  operation  strategies  for 
maximizing  the  plant  revenue  should  be  defined. 

In  such  a  context,  several  studies  have  been  carried  out  to 
define  the  optimal  bidding  strategy  for  individual  PHES  with 
reference  to  different  time  horizons  (daily,  weekly  or  monthly). 
All  these  algorithms  were  based  on  an  objective  function  to  be 
maximized,  taking  into  account  some  technical  constraints. 

In  2004,  Lu  et  al.  [38]  focused  their  attention  on  an  individual 
PHES  with  a  storage  capacity  of  more  than  10  h  and  developed  a 
multi-stage  optimization  algorithm  based  on  a  weekly  horizon, 
considered  to  be  a  proper  period  for  the  optimal  use  of  the  water 
storage  capacity.  Their  algorithm  maximizes  the  profit  of  the  plant 
by  considering  the  revenue  deriving  by  both  the  day-ahead  market 
and  the  ancillary  services  market.  According  to  their  model,  the  day- 
ahead  market  price  is  forecasted  on  an  hourly  basis  in  a  weekly  time 
period,  whereas  the  ancillary  market  prices  (Pre  and  Pm)  are  supposed 
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to  be  constant 

tg  _  fp 

k[0,T]  —  2  EgenPgen(i)4-ErsPrstrs4-ErnPrntrn  ~  X  EpumpPpump(t)  Co&M  (2) 

i  =  1  1=1 

In  the  optimization  algorithm,  the  weekly  forecasted  day-ahead 
market  curve  (Fig.  5a)  is  sorted  by  ascending  prices  (Fig.  5b)  and 
the  profit  of  the  PHES  plant  is  optimized  by  progressively  increasing 
the  pumping  time  tp  and  checking  for  the  maximum  profit  (Eq.  (2)). 

The  technical  constraint  represented  by  the  maximum  storable 
energy  due  to  the  upper  reservoir  dimensions  (£max)  is  not 
considered  in  the  optimization  algorithm,  but  it  is  applied  as  a 
corrective  criterion  of  the  optimal  solution.  In  the  limited  time 
intervals  of  the  optimal  solution  in  which  the  constrain  is  violated, 
the  unconstrained  optimization  algorithm  is  performed  again  until 
no  violations  are  identifiable  in  the  weekly  time  period  (Fig.  6). 

In  2009,  Kazempour  et  al.  [20]  proposed  a  weekly-based 
optimization  algorithm  still  aimed  at  the  maximization  of  the 
plant  profit.  This  algorithm  considers  the  weekly  forecast  not  only 
of  the  daily-ahead  market  but  also  of  the  ancillary  markets  and  in 
particular  of  the  spinning  and  regulation  markets  on  an  hourly- 
based  price  variation  (Eq.  (2)). 

In  2011,  Connolly  et  al.  40]  limited  their  analysis  to  the  day- 
ahead  market  and  compared  different  operation  strategies  aimed 
at  maximizing  the  plant  revenue  deriving  from  the  electricity 
trading  in  one  year  of  plant  operation  (T=  8760  h) 

tg  _  fp 

k[0,q  =  2  EgcnPgen(i)  "  ^  EpumpPpumpO)  ~  Co&M  (3) 

i = 1  1=1 


€120 


Fig.  6.  Summarizing  scheme  of  the  optimization  strategy  proposed  by  Lu  et  al.  [38], 


They  demonstrated  that  in  the  hypothesis  of  an  accurate  predic¬ 
tion  of  the  day-ahead  hourly  prices,  an  operating  strategy  aimed  at 
maximizing  the  revenue  only  of  the  next-day  can  be  more  efficient 
than  an  operation  strategy  aimed  at  maximizing  the  revenue  of 
the  entire  year. 

Algorithms  were  also  proposed  to  define  optimal  operation 
strategies  of  combined  hybrid  wind/photovoltaic  hydro  energy 
storage  plants.  In  2004,  Castronuovo  et  al.  [39]  developed  an 
optimization  approach  by  means  of  a  primal-dual  interior  point 
algorithm  aimed  at  maximizing  the  plant  profit  £[0,r]  on  a  2-day 
horizon  period  (T=48  h) 

tg  _  h 

k[0T|  =  X  EgenPgen0)~  X  EpumpPpump(t)  ~  Cs[AEcap  -  -  Cp  Afco  (4) 
i = 1  (=1 


where  Cst  and  CE  are  properly  defined  penalties  for  increases  in 
storage  capacity  A Ecap  and  in  the  initial  energy  level  A£0  of  the 
upper  reservoir. 

In  this  case,  the  plant  production  Egen  is  the  sum  of  the  energy 
produced  by  the  hydro  generator£flJ,(iro|genand  of  the  amount  of 
energy  produced  by  the  wind  generator  and  delivered  to  the 
network  Ewind\gen 


Egen  —  Ewind  \„en  +  £  hydro  L 


(5) 


b 


Time  [hr] 

Fig.  5.  Convenience  limits  for  the  operating  modes  of  a  PHES  in  the  day-ahead 
market:  (a)  weekly  day-ahead  market  curve;  (b)  weekly  day-ahead  market  curve 
sorted  by  ascending  prices. 


whereas  Epump  is  the  energy  purchased  from  the  network  for 
pumping  water  that,  added  to  the  possible  surplus  of  wind 
generation  not  delivered  to  the  network  Ewind  |  ,  defines  the 

global  amount  of  energy  stored 

Estorage  =  £wind \pump  +  Epump  (6) 

As  regards  the  reservoir  capacity,  the  model  makes  it  possible  to 
evaluate  the  minimum  amount  of  incremental  storage  capacity  to 
follow  a  specified  power  demand  curve  and  thereby  also  optimize 
the  storage  size.  As  regards  technical  constraints,  a  maximum 
allowed  wind  power  exchange  with  the  system  was  fixed  to 
guarantee  the  grid  stability.  Even  though  the  strategy  definition 
only  considers  the  day-ahead  market,  the  combined  optimization 
of  strategy  and  storage  sizing  made  it  possible  to  gain  not  only  in 
plant  profit,  but  also  in  wind  energy  production.  The  importance  of 
an  optimal  sizing  of  a  combined  wind-hydro  storage  plant  was 
studied  in  depth  by  Khatibi  and  Jazaeri  [41],  that  presented 
different  optimization  techniques  (single  or  multi-objective 
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optimizations,  parametric  studies,  sensitivity  studies)  aimed  at 
maximizing  the  wind  energy  penetration  by  optimal  plant  design. 

The  possibility  of  increasing  the  wind  energy  penetration 
through  the  optimal  combination  with  PHES  was  also  demon¬ 
strated  by  Dursun  et  al.  [42]  and  by  Bueno  and  Carta  [43]  that 
developed  a  techno-economical  multi-stage  algorithm.  On  the 
basis  of  technical  and  economic  databases  of  the  commercial  plant 
components,  this  algorithm  defines  an  optimum-sized  system 
aimed  at  maximizing  the  exploitation  of  the  wind  energy  source 
without  penalizing  the  grid  stability.  The  model  considers  not 
only  investment,  operating  and  maintenance  costs,  but  also  those 
costs  associated  with  environmental  damage  caused  by  energy 
production. 

In  2012,  Ding  et  al.  [44]  presented  a  stochastic  optimization  of 
the  daily  operation  strategy  of  a  hybrid  wind-hydro  energy  plant 
aimed  at  maximizing  plant  revenue.  To  consider  the  influence  of 
forecast  uncertainties  of  the  next-day  wind  power  output  on  the 
optimization,  the  model  was  based  on  chance  constraints  and 
scenario  analyses.  Start-up  and  shutdown  costs  of  the  PHSP  units 
as  well  as  penalties  for  output  deviation  were  considered. 

Much  more  complex  algorithms  have  been  recently  proposed 
to  formulate  an  efficient  management  optimization  procedure  of 
large-scale  hydropower  systems,  aimed  at  increasing  the  exploita¬ 
tion  of  the  intermittent  renewable  sources  and  at  supporting  and 
increasing  the  power  grid  stability.  However,  the  operation  of 
these  systems  represents  a  mixed,  non-linear  and  non-convex 
programming  problem  and,  up  till  now,  optimal  solutions  have 
only  been  found  by  approximation  of  the  original  problem. 
Different  methods  were  adopted:  dynamic  programming  models 
[45],  non-linear  programming  models  [46],  genetic  algorithms 
[47],  simulated  annealing  [48],  ant  colony  algorithms  [49],  multi¬ 
stage  combined  models  [50],  and  particle  swarm  optimization 
models  [51-55], 

These  algorithms,  in  comparison  with  the  above  mentioned 
simplified  ones,  have  certainly  increased  the  possibility  of  con¬ 
sidering  a  complex  model  by  increasing  the  number  of  variables  to 
optimize.  However,  the  definition  of  a  model  properly  considering 
the  techno-economical  interaction  between  the  electricity  market 
and  the  operating  plant  capability  still  remains  a  key-challenge  not 
only  for  working  plants  in  terms  of  revenue  maximization  but  also 
for  planned  plants  since  it  allows  to  evaluate  the  investment 
profitability  on  one  side,  and  to  define  the  most  proper  design 
strategy  on  the  basis  of  the  electricity  market  in  which  the  plant 
will  be  required  to  operate  on  the  other  side  (as  seen  in  Section 
2.1).  For  example,  as  regards  the  electricity  market,  to  consider 
only  the  day-ahead  market  or  a  simplified  model  of  the  ancillary 
markets  could  negatively  affect  the  evaluation  of  the  investment 
profitability  of  a  PHES  since  in  the  next  10  years  the  forecasted 
spread  between  the  maximum  and  minimum  electricity  price  in 
the  day-ahead  market  will  be  significantly  reduced  due  to  the 
flattening  of  daily  power  load  curve.  As  a  consequence  of  the 
flattening  of  the  electricity  price  curve,  in  a  long-term  view  the 
revenue  for  a  PHES  will  mainly  derive  from  the  ancillary  services 
and  the  correct  evaluation  of  the  plant  profitability  cannot  leave 
these  revenues  out  of  consideration  or  simplify  their  determina¬ 
tion.  Another  important  aspect,  frequently  neglected  in  the  mod¬ 
els,  is  the  influence  of  the  forecasts  uncertainties  of  the  electricity 
markets  and  of  the  wind/PV  power  output  on  the  analysis  results. 

To  correctly  evaluate  the  capability  of  a  PHES  of  interacting 
with  the  electricity  markets,  the  economical  model  of  the  markets 
should  be  properly  combined  with  an  operational  model  of  the 
plant  that  considers  the  real  grid-connection  capability  of  the 
plant  by  taking  into  account  head  and  flow  rate  operating  range, 
variable-speed  effects  on  the  operating  range,  pump-turbine 
stability  limits,  storage  capacity,  losses,  response  times  for  load 
variation  in  generation  and  pumping  mode,  response  times  for 


variation  of  the  working  mode  (pumping  to  production  and 
production  to  pumping)  and  dynamical  aspects. 

Up  till  now,  none  of  the  studies  proposed  in  literature  have 
considered  in  details  all  the  above  mentioned  economical  and 
technical  aspects  and  hence  the  ambitious  goal  of  a  realistic 
simulation  of  the  interaction  between  the  PHES  and  the  electricity 
market  still  remains  a  challenge  due  to  the  complexity  of  system 
to  model  and  due  to  the  significant  computing  costs. 


3.  Small  hydro  power 

Hydropower  on  a  small  scale  has  not  an  international  agreed 
definition  and  its  upper  limit  can  vary  between  2.5  and  30  MW, 
even  though  the  most  widely  accepted  value  is  10  MW. 

Being  among  the  energy  storage  technologies  currently  avail¬ 
able  on  a  small  scale,  it  is  recognized  as  one  of  the  most  cost- 
effective  because  of  its  predictable  energy  characteristics,  its  long 
term  reliability  and  its  reduced  environmental  effects  [56,57], 
Unlike  large  scale  hydropower,  whose  development  is  associated 
with  the  building  of  large  dams,  small  hydro  power  plants  is 
generally  “run-of-river”,  with  a  weir  crossing  the  river  and  little  or 
no  water  stored. 

Like  pumped-hydro  energy  storage  plants,  small  hydro-power 
plants  require  specific  characteristics  of  the  installation  sites  to  be 
feasible,  such  as  a  natural  head  and  the  possibility  of  water  power 
harnessing  [58].  Depending  on  site  characteristics,  such  as  avail¬ 
able  head  and  flow  rate,  and  on  the  selected  running  speed  of  the 
generator,  different  turbine  types  with  different  operating  range 
and  performance  characteristics  are  available  and  can  be  divided 
into  impulse  and  reaction  turbines. 

In  impulse  turbines  the  water  remains  at  atmospheric  pressure 
and  the  energy  exchange  is  only  based  on  the  kinetic  energy 
variation,  whereas  reaction  turbines  are  characterized  by  a  runner 
fully  immersed  in  water.  Since  reaction  turbines  exploit  the 
available  head  in  terms  of  both  pressure  and  kinetic  energy 
variation,  a  draft  tube  is  placed  downstream  the  runner  not  only 
to  discharge  the  fluid  but  also  to  increase  the  energy  exchange 
owing  to  a  further  pressure  reduction  induced  at  the  runner  outlet. 

Impulse  turbines  can  be  divided  in  Pelton,  Turgo  and  cross-flow 
turbines.  The  most  important  impulse  turbine  is  undoubtedly  the 
Pelton  turbine,  patented  in  1880  by  Lester  Allen  Pelton.  It  is 
composed  by  one  or  more  penstocks,  equipped  with  proper 
nozzles,  and  by  a  wheel,  whose  blades  are  split  buckets  properly 
distributed  around  its  periphery  (Fig.  7a).  The  water  jet,  coming 
out  from  each  nozzle  and  directed  tangentially  at  the  wheel,  heats 
the  buckets  and  is  split  in  half,  before  being  deflected  back  almost 
180°.  The  consequent  moment  of  momentum  variation  determines 
the  turbine  energy  exchange. 

Turgo  turbines  are  similar  to  Pelton  turbines,  but  are  character¬ 
ized  by  a  different  blade  shape  and  by  a  different  relative 
dispositions  between  runner  and  penstocks  (Fig.  7b).  The  jet 
strikes  the  plane  of  the  runner  at  an  angle  of  about  20°  so  that 
the  water  enters  the  runner  from  one  side  and  comes  out  from  the 
other,  eliminating  the  interaction  hazard  between  flows  dis¬ 
charged  by  different  buckets,  a  risk  that  could  affect  the  Pelton 
operating  range  instead,  thus  limiting  the  maximum  flow  rate. 

Cross-flow  turbines  are  simpler  and  cheaper  than  the  other 
impulse  turbines  and  were  patented  after  some  evolutions  in  1933 
by  A.S.  Mitchell  and  F.  Ossberger  (Fig.  7c).  The  jet  water  enters  on 
one  side  of  a  drum-like  runner  and  leaves  it  on  the  opposite  side 
after  a  double  passage  through  the  runner  blades. 

Reaction  turbines  can  be  divided  in  Francis  and  axial  turbines. 
Francis  turbines,  patented  by  James  Bickens  Francis  at  the  end  of 
the  nineteenth  century,  are  mixed-flow  reaction  turbines  com¬ 
posed  by  spiral  case,  distributor,  runner  and  draft  tube  (Fig.  7d). 
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Fig.  7.  Scheme  of  different  types  of  hydraulic  turbines:  (a)  single-jet  Pelton  turbine;  (b)  Turgo  turbine;  (c)  cross-flow  turbine;  (d)  Francis  turbine;  and  (e)  Kaplan  turbine. 


The  spiral  case  receives  the  flow  and  distributes  it  peripherally  to 
the  distributor,  represented  by  mobile  guide  vanes  whose  aim  is  to 
control  the  flow  rate  entering  into  the  runner  and  to  define  the 
flow  angle  at  the  runner  inlet. 

The  axial  flow  reaction  turbines  include  Kaplan  and  propeller 
turbines  and  are  generally  composed  by  guide  vanes,  a  runner  and 
a  draft  tube.  Once  it  entered  the  turbine,  the  water  received  a  swirl 
component  by  the  guide  vanes,  that  is  almost  completely  absorbed 
by  the  runner,  thus  determining  the  energy  exchange.  The  water  is 
often  directed  in  the  draft  tube  with  a  little  residual  angular 
momentum  to  prevent  boundary  layer  detachments.  Propeller 
turbines  are  often  characterized  by  adjustable  guide  vanes  to  vary 
the  flow  rate,  whereas  in  Kaplan  turbines  even  the  runner  blade 
can  be  adjusted  (Fig.  7e).  Even  if  the  mechanisms  for  adjusting 
runner  and  guide  vanes  increase  the  machine  cost,  they  bring 
about  a  significant  improvement  of  the  efficiency  at  part  load,  as 
can  be  seen  in  Fig.  8. 

The  comparison  in  terms  of  efficiency  and  operating  range 
between  turbines  highlights  that  Pelton,  Kaplan  and  cross-flow 
turbines  retain  high  efficiency  values  even  at  part  loads;  instead, 
quick  Francis  turbines  and  fixed  pitch  propeller  turbines  are 
characterized  by  a  rapid  efficiency  fall  below  50%  and  80%  of  the 
design  flow  rate  respectively. 


3.1.  Optimal  sizing 

Since  most  of  the  rare  high-head  installation  sites  have  been 
already  exploited,  the  great  potential  of  small  hydro  power  is 
represented  by  low-head  sites,  whose  economical  attractiveness  is 
strictly  connected  with  an  optimal  design  of  the  electro-mechanical 
equipment,  mainly  represented  by  turbine  and  generator. 


The  capacity  sizing  of  a  small  hydropower  plant  is  strictly 
connected  with  flow  availability  and  is  based  on  the  analysis  of  the 
flow  duration  curve  (Fig.  9).  The  flow  duration  curve  is  a  cumu¬ 
lative  distribution  of  the  stream  flow  rate  of  a  site  on  an  annual 
basis  and  shows  the  percentage  of  time  that  the  site  water  flow 
equals  or  exceeds  a  specific  value. 

The  available  flow  rate  range  of  the  site  is  included  between  the 
maximum  flow  rate  flowing  for  at  least  1  day  per  year  (Qmax)  and 
the  minimum  flow  rate  flowing  into  the  site  (Qmm) 

Q-min  <  Q  <  Qmax  O) 

This  flow  rate  cannot  be  fully  exploited,  since  it  is  established  by 
law  that  a  reserved  flow  (RF),  a  sort  of  environmental  flow,  should 
be  released  downstream  to  keep  the  ecosystem  in  the  conditions 
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Fig.  9.  Example  of  flow  duration  curve. 


that  prevaled  before  the  power  plant  construction.  Therefore,  the 
flow  rate  range  effectively  exploitable  is  obtained  from  the  flow 
duration  curve  after  the  RF  subtraction. 

A  further  limit  on  the  flow  rate  exploitation  is  represented  by 
the  turbine  chosen  since  each  turbine  operates  with  acceptable 
efficiency  between  a  minimum  flow  rate  QS,m  and  a  maximum 
flow  rate  QJj,ax  (Fig.  9).  Therefore,  for  flow  rates  smaller  than  QJ!n 
the  plant  should  be  stopped,  releasing  all  the  available  flow  rates, 
whereas  for  flow  rates  greater  then  Q,*  ax  the  plant  is  not  be  able  to 
exploit  the  exceeding  water  quantity  released  downstream  with¬ 
out  being  exploited.  Depending  on  the  turbine,  the  energy  pro¬ 
duction  over  a  time  period  T  can  be  determined  as 

E=  lo  P{t)dt=  L  pgQ{t)h{t)’lhydit)’lmechrie>dt  (8) 

where  P(t)  is  the  power  obtained  by  the  turbine,  p  is  the  water 
density,  Q(t )  is  the  flow  rate,  /i(t)  is  the  net  head,  tjhyd  is  the 
hydraulic  efficiency  (Fig.  8),  i}meCh  and  //ei  are  respectively  the 
mechanical  and  the  electrical  efficiencies.  The  net  (also  called 
“effective”  or  “available”)  head  h(t)  is  obtained  by  the  gross  head 
hg  (i.e.  the  difference  between  the  level  of  the  upper  reservoir  and 
the  level  of  the  tailwater)  after  subtraction  of  the  hydraulic  losses 
hft)  in  the  pipeline  connecting  the  reservoir  with  the  turbine  and 
of  the  kinetic  energy  discharged  by  the  turbine  without  being 
exploited 

h(t)  =  hg-(hr(t)+Vk^J  (9) 

where  vouClet  is  the  water  velocity  at  the  draft  tube  outlet. 
According  to  Eq.  (9),  the  draft  tube  is  an  integral  part  of  the 
turbine. 

Several  studies  focused  their  attention  on  the  design  of  a  small 
hydro  power  plant  as  its  features  played  a  key  role  in  the  definition 
of  the  plant  profitability.  Flosseini  et  al.  [59]  determined  the 
optimal  installation  capacity  of  a  plant  as  a  compromise  between 
technical,  economic  and  reliability  indices.  However,  the  algorithm 
did  not  enter  into  the  details  of  the  electro-mechanical  design. 

A  non-linear  programming  method  was  adopted  by  Almeida  et  al. 
[60]  who  optimized  the  project  configuration  of  a  small-hydro  power 
plant  by  considering  technical,  hydrological,  operational,  economic 
and  financial  aspects.  The  model  considered  different  possible 
hydrological  and  electricity  market  scenarios,  that  significantly 
affected  the  optimal  project  configuration.  In  2007,  Anagnostopoulos 
and  Papantonis  [61  ]  carried  out  a  multi-objective  optimization  study 
to  determine  the  best  compromise  between  the  economic  benefits  of 
the  investment  and  the  hydraulic  potential  exploitation  in  terms  of 


annual  energy  production  and  river  water  exploitation.  Niadas  and 
Mentzelopoulos  [62]  focused  their  attention  on  the  influence  of  the 
flow  duration  curve  on  plant  capacity  sizing  and  presented  a 
methodology  for  its  probabilistic  treatment.  The  influence  of  the  flow 
duration  curve  on  plant  sizing  was  also  highlighted  by  Pena  et  al.  [63] 
who  developed  a  procedure  for  estimating  the  water  flow  of  a  site  on 
the  basis  of  time  series  forecasting  methods.  Different  techniques 
aimed  at  enhancing  the  available  historical  data  were  verified  and  the 
obtained  flow  duration  curve  was  exploited  for  determining  the 
turbine  design  flow  rate  on  the  basis  of  a  fixed  capacity  factor.  A 
different  approach  was  proposed  by  Borges  and  Pinto  in  2008  [64], 
who  presented  a  stochastic  method  for  analysing  the  influence  of 
river  inflow  variations  and  generator  unit  availability  on  the  power 
plant  availability.  Their  analysis,  starting  from  known  plant  character¬ 
istics  (historical  river  flow  data,  height,  turbine  type,  turbine  char¬ 
acteristics,  generator  characteristics),  makes  it  possible  to  determine 
the  expected  value  of  the  annual  power  plant  generation,  the  flow 
duration  curve  and  several  reliability  indices.  Possible  generation  unit 
failures  are  also  considered  in  order  to  define  also  the  risk  level 
associated  with  the  expected  energy  production.  In  2008,  Bockman 
et  al.  [65]  combined  the  optimal  sizing  of  a  small  hydropower  plant 
with  an  analysis  of  the  influence  of  the  uncertain  electricity  prices  on 
the  plant  revenues  and  developed  a  real  options-based  model  to  find 
the  electricity  price  limit  for  a  profitable  investment. 

In  2011,  Santolin  et  al.  [66]  considered  the  importance  of 
significant  technical  aspects,  such  as  turbine  type,  cavitation 
behaviour  and  turbine  dimensions,  in  determining  the  capacity 
installation  size  and  proposed  a  method  based  on  technical  and 
economic  aspects.  Seven  design  parameters  were  considered:  the 
turbine  type,  the  main  turbine  dimensions,  the  maximum  installa¬ 
tion  height  to  avoid  cavitation  inception,  the  annual  energy 
production,  the  installation  costs,  the  Net  Present  Value  (NPV) 
and  the  Internal  Rate  of  Return  (IRR)  as  economic  indices.  Starting 
from  the  flow  duration  curve,  the  procedure  highlighted  the 
importance  of  a  simultaneous  analysis  of  all  these  aspects  to  a 
proper  definition  of  the  design  operating  conditions  aimed  at 
guaranteeing  not  only  plant  profitability,  but  also  its  feasibility. 

Due  to  the  last  advances  in  hydro  technology,  one  important 
aspect  to  consider  in  the  analysis  of  the  optimal  design  in  small 
hydro  applications  should  be  the  use  of  the  variable  speed 
technology  in  installation  sites  characterized  by  significant  head 
and/or  flow  rate  variations.  The  possible  application  of  this 
technology  could  significantly  increase  the  plant  profitability  by 
increasing  the  hours  of  production,  and  its  implementation  in  the 
optimization  algorithm  will  certainly  represent  a  significant 
advance  in  the  optimal  sizing  of  a  small  hydro  power  plant. 
Moreover,  since  the  uncertainties  in  the  flow  duration  curve  and 
in  the  forecast  electricity  prices  could  affect  the  analysis  results,  a 
risk  analysis  related  to  these  aspect  could  represent  a  future 
challenge  to  add  to  standard  optimization  strategy. 

3.2.  Optimal  operation  strategy  for  small  hydro  power 

The  choice  of  proper  design  operating  conditions  is  undoubt¬ 
edly  a  fundamental  goal  for  maximizing  both  investment  profit¬ 
ability  and  the  exploitation  of  the  hydraulic  potential  of  water 
resources.  However,  once  the  plant  goes  into  operation,  to  main¬ 
tain  the  unit  efficiency  and  to  maximize  the  plant  revenue,  the 
electro-mechanical  component  performance  should  be  preserved 
and  the  operating  strategy  should  be  optimized  in  relation  to 
inflow  variations  and  market  fluctuations. 

In  2003,  Liu  et  al.  [67]  developed  a  method  of  monitoring  and 
verifying  the  economic  performance  of  an  operating  hydropower 
plant  by  using  new  performance  indices  based  on  the  reachable 
efficiency  and  the  operating  efficiency  associated  with  the  main¬ 
tenance  level  and  with  the  operation  management  level  of  the 
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plant  respectively.  These  new  concepts  provide  decision  making 
support  tools,  both  for  planning  refurbishment  and  maintenance 
activities  and  for  suggesting  new  operation  strategies. 

To  analyse  the  plant  performance  under  specified  operating 
conditions,  Garrido  et  al.  [68]  developed  a  simulation  tool  based 
on  an  object  oriented  modelling  language,  that  allowed  to  model 
different  configurations  of  hydropower  stations  with  the  possibi¬ 
lity  of  choosing  different  numbers  of  turbines  and  spillway  gates. 
The  tool  made  it  possible  configure  plant  parameters  and  to 
simulate  the  plant  behaviour  at  different  operating  conditions 
and  thereby  obtain  useful  information  on  the  plant  behaviour, 
such  as  reservoir  level,  water  flows,  turbine  efficiency  and  so  on. 
However,  simulation  tools  like  those  proposed  by  Liu  et  al.  [67]  are 
highly  user  dependent  since  they  are  based  on  the  plant  para¬ 
meters  manually  adjusted  by  the  operator  and  do  not  guarantee 
the  definition  of  an  optimal  operating  strategy,  such  as  that 
obtained  by  applying  a  properly  defined  optimization  procedure. 

Fleten  and  Kristoffersen  [69]  developed  a  multi-stage  mixed- 
integer  linear  stochastic  program  to  define  an  optimized  short-term 
production  plan  complying  with  the  day-ahead  commitments  of  the 
previous  day.  Starting  from  a  production  plan  scheduled  in  time  by 
the  results  of  the  day-ahead  auction,  the  procedure  sets  out  to 
identify  the  expected  production  between  the  various  parts  of  the 
plant  in  order  to  achieve  effective  and  efficient  operation.  The 
uncertainties  of  the  day-ahead  market  prices  and  of  the  reservoir 
inflows  due  to  the  variable  market  conditions  and  to  the  unpredict¬ 
able  weather  situations  are  taken  into  consideration  in  the  planning 
by  the  stochastic  programming  framework.  The  procedure  was 
originally  proposed  for  small  hydro  power  plants,  but  could  be 
implemented  in  more  complex  large  hydropower  systems 
thanks  to  the  mixed-integer  linear  programming  approach  that 
makes  it  possible  to  model  start-up  costs  and  discrete  hydro  unit- 
commitment  constraints. 

A  further  improvement  in  the  definition  of  an  optimized  short¬ 
term  production  plan  was  proposed  by  Catalao  et  al.  [70],  whose 
evaluation  was  based  on  a  mixed-integer  non-linear  programming 
approach  considering  the  non-linear  dependence  between  power 
generation,  water  discharge  and  head.  Even  though  this  approach 
was  shown  to  increase  the  computational  costs  of  the  analysis,  it 
was  possible  to  define  hydroelectric  operating  conditions  with 
greater  accuracy  and,  in  particular,  to  take  into  account  the  head 
dependency  of  a  short-term  production  plan. 

Even  though  these  algorithms  represent  a  significant  advance 
in  the  definition  of  the  optimal  management  strategies  of  a  small 
hydro  power  plant,  they  require  to  be  further  improved  by  the 
implementation  of  models  considering  series/parallels  of  large/ 
small  hydro  power  plant  in  which  the  management  strategy  of  one 
plant  is  affected  by  the  management  of  the  others.  An  example  is 
represented  by  one  of  the  most  interesting  possible  application  of 
the  variable-speed  technology  in  the  small  hydro  field,  that  is  the 
exploitation  of  the  reserved  flow  of  a  large  hydro  power  plant  [71  ]. 
The  modelling  of  the  head  dependency  of  the  small  hydro 
production  plan  from  the  management  strategy  of  the  large  one 
(that  could  be  a  PHES)  is  undoubtedly  an  extremely  complex  aim 
and  is  one  of  the  most  interesting  challenge  in  this  research  field. 


4.  CFD  as  a  tool  for  improving  hydraulic  turbine  and  pump- 
turbine  performance 

The  peculiarity  of  Computational  Fluid  Dynamics  (CFD)  is  its 
capacity  to  calculate  the  flow  field  around  an  arbitrary  obstacle  or 
through  a  channel  of  an  arbitrary  shape  via  the  numerical  solution 
of  the  equations  that  describe  the  flow  development  (Navier- 
Stokes  equations,  energy  equation,  global  and  partial  continuity 
equations,  etc.). 


Over  the  years,  due  to  the  limits  in  calculation  capacity  and 
speed,  not  every  type  of  flow  field  (unsteady,  three  dimensional, 
compressible,  turbulent,  reacting,  etc.)  could  be  calculated  and  the 
governing  equations  were  simplified  in  order  to  describe  with 
more  or  less  accuracy  the  flow  field  to  be  analysed. 

The  history  of  computational  fluid  dynamics  started  at  the 
beginning  of  the  last  century  when  the  reduced  (first  human  and 
then  numerical)  computation  capacity  limited  the  analyses  on 
two-dimensional  flow  fields,  whose  governing  equations  were 
simplified  by  the  linearized  potential  equations  solved  by  con¬ 
formal  mapping  or  singularity  method  [72,73],  It  was  only  in  the 
1970s  that  the  first  3D  models  were  numerically  analysed  by 
application,  in  the  fluid-dynamics  field,  of  the  finite  difference 
method  (FDM)  and  the  finite  element  method  (FEM),  historically 
originating  from  the  field  of  structural  mechanics  and  aerody¬ 
namics  [74,75],  These  codes  demonstrated  the  possibility  of 
numerically  solving  with  good  accuracy  both  2D  and  3D  potential 
flow  problems  and  the  possibility  of  analysing  complex  geometries 
thanks  to  the  FEM  way  of  modelling  by  triangular  and  rectangular 
8-node  parabolic  elements. 

As  regards  FEM  applications  in  the  turbomachinery  field,  fluid- 
dynamical  analyses  of  stationary  and  rotating  rows  were  initially 
carried  out  by  means  of  quasi-3D  approaches  in  which  the  flow 
was  approximately  reconstructed  by  a  combination  of  the  numer¬ 
ical  meridional  and  blade-to-blade  flow  fields  [76,77].  In  spite  of 
the  limitations  of  such  an  approach,  whose  results  were  reason¬ 
ably  accurate  only  for  strongly  accelerated  flows  and  close  to  the 
best  efficiency  point,  some  useful  information  could  be  already 
derived  for  turbomachine  designers. 

However,  even  though  improvements  to  the  quasi-3D  inviscid 
approach  were  successively  obtained  by  adding  3D-boundary  layer 
calculations,  the  simulation  of  the  flow  fields  in  radial  runners  was 
still  not  accurate  enough  since  the  majority  of  the  3D-effects 
developing  in  them  was  driven  by  vorticity  rather  than  by 
viscosity.  For  these  reasons,  in  the  1990s,  the  need  to  correctly 
model  the  vorticity-driven  secondary  flows  resolutely  addressed 
towards  a  fully  inviscid  3D  Euler  code,  as  the  computational 
capacity  and  the  method  accuracy  were  not  still  ripe  for  directly 
stepping  into  Navier-Stokes  simulations  [72].  These  codes  made  it 
possible  to  analyse  with  a  high  degree  of  accuracy  the  runner 
performance  at  design  and  off-design  operating  conditions,  as  long 
as  the  real  flow  was  not  dominated  by  viscous  effects,  and  this 
capability  made  it  possible  to  infer  the  real  possibility  of  trans¬ 
forming  the  numerical  analysis  from  a  research  tool  to  a  practical 
and  extremely  useful  design  tool.  However,  in  this  period  the 
hydraulic  industry  was  facing  important  challenges  due  to  the 
emergence  of  a  huge  demand  for  rehabilitation  and  upgrading  of 
old  hydropower  plants.  Relying  only  on  previous  experience  by 
scaling  previous  turbine  designs  was  becoming  a  hazardous 
strategy  owing  to  the  numerous  constrains  imposed  by  the 
existing  plants  and  the  uncertainties  connected  to  the  attainment 
of  the  increase  in  demand  of  power.  Furthermore,  resorting  to 
model  tests  was  also  too  expensive  because  the  cost  of  a  medium 
size  runner  model  was  nearly  the  same  as  the  full-scale  prototype. 
So,  to  face  the  increased  demand  for  very  high-performance 
hydropower  plants,  new  CFD-based  design  approaches  were 
proposed,  in  which  the  new  runners  were  supposed  to  be  analysed 
and  compared  with  the  old  ones  by  CFD,  greatly  limiting  the 
support  of  model  tests.  This  revolutionary  design  strategy  pro¬ 
moted  the  introduction  of  Navier-Stokes  codes  for  the  computa¬ 
tion  of  both  turbulent  and  viscous  effects.  In  contrast  to  Euler 
equations,  Navier-Stokes  equations  allowed  the  detection  of  the 
boundary  layer  separation  in  all  turbine  components  as  well  as  the 
loss  analysis  both  at  design  and  off-design  operating  conditions. 
By  2000,  the  solution  of  the  steady  Reynolds  averaged  Navier- 
Stokes  (RANS)  equations,  based  on  the  finite  volume  modelling 
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methods,  was  demonstrated  to  be  the  only  tool  for  realistically 
analysing  the  onset  of  pump-turbine  instabilities  as  well  as  of 
viscosity-based  secondary  flows  and  hence  became  the  state-of- 
art  in  the  major  refurbishment  projects  [78]. 


4.1.  Interface  model  between  rotating  and  stationary  components 

The  comparison  between  numerical  and  experimental  results 
are  strongly  influenced  by  boundary  conditions,  which  are  not  at 
all  easy  to  set  when  different  components,  rotating  and  stationary, 
are  fluid-dynamically  coupled  to  each  other. 

Test  data  of  a  Francis  turbine  [79]  have  demonstrated  that  the 
pronounced  low  pressure  peak  at  runner  inlet  cannot  be  revealed 
by  CFD  analyses  when,  as  in  the  case  of  boundary  conditions  at  the 
runner  inlet,  the  experimental  inlet  velocity  profiles  (circumfer¬ 
entially  averaged)  are  adopted.  Instead,  a  good  adherence  between 
numerical  and  experimental  tests  may  be  observed  when  a 
stacking  technique  between  stationary  upstream  guide  vanes 
and  the  rotor  is  considered  (Fig.  10)  and  nowadays  this  technique 
represents  a  well-consolidated  approach  for  numerical  studies  on 
cavitating  flows  in  Francis  turbines  [80], 

During  the  mid-80s  the  development  of  interface  models 
between  rotating  and  stationary  domains  made  it  possible  to 
appraise  the  interaction  between  adjacent  turbomachinery  com¬ 
ponents,  opening  the  era  of  CFD  simulations  of  entire  turbines. 
This  interaction  was  modelled  by  more  or  less  sophisticated 
mathematical  algorithms,  whose  common  peculiarity  was  to 
reduce  the  time-dependent  interaction  between  rotating  and 
stationary  components  to  a  steady-state  problem.  These  algo¬ 
rithms  modelled  in  a  different  way  the  transportation  of  the 
stator/rotor  blade  wakes  across  the  interface,  resulting  in  an 
interface  represented  by  a  circumferential  mixing  plane  or  by  a 
locally  non-uniform  flow  field,  depending  on  the  algorithm. 
However,  unlike  the  transient  interface  approach,  blade  wakes 
driven  by  the  flow  through  a  stator-rotor  interface  could  not  be 
detected  adequately  when  a  steady  frozen-rotor  interface  algo¬ 
rithm  was  adopted,  owing  to  the  possible  occurrence  of  over¬ 
estimations  of  the  wake  disturbance  effects  [81], 

The  ability  and  capability  of  modelling  the  flow  through  the 
entire  machine  in  a  single  CFD  simulation  further  boosted  the  use 
of  the  CFD  analysis,  favouring  its  evolution  from  a  design  tool  into 
a  high-level  analysis  tool,  as  it  is  today,  to  be  used  in  all  stages  of 
the  refurbishment,  from  the  feasibility  study  to  the  design  and 
optimization  phase. 


In  spite  of  these  significant  advances,  the  need  of  investigating 
the  unstable  behaviour  of  turbine  and  pump-turbine  to  enlarge 
their  operating  range  will  require  further  improvements  in  CFD  in 
order  to  model  turbulent  flow  and  performance  at  a  wider  range  of 
discharge  rates  more  accurately  and  faster. 

4.2.  Pelton  CFD 

In  recent  years,  significant  developments  in  CFD  testified  by  the 
accurate  results  of  the  numerical  analysis  on  Francis  and  Kaplan 
turbines,  and  paved  the  way  for  the  exploitation  of  CFD  even  to  the 
more  complex  Pelton  turbines. 

Unlike  the  Francis  and  Kaplan  turbines,  whose  design  procedures 
are  well-consolidated  also  thanks  to  CFD,  Pelton  turbines  are  still 
characterized  by  semi-empirical  design  criteria  because  of  the  diffi¬ 
culty  in  investigating  (both  experimentally  and  numerically)  the 
complex  fluid-dynamical  interaction  between  the  water  jet  and  the 
rotating  buckets,  which  greatly  affects  turbine  performance. 

With  the  last  improvements  in  computational  accuracy  and 
efficiency,  the  experimental  investigations  were  successfully  com¬ 
bined  with  numerical  analyses  to  gain  new  insights  in  such 
complex  phenomena. 

Following  some  attempts  at  studying  a  moving  Pelton  turbine 
bucket,  it  was  only  in  2006  that  Perrig  et  al.  [82]  succeeded  in 
performing  a  numerical  analysis  of  the  cutting  process  of  axial- 
symmetric  jets  on  the  bucket  cut-out,  whereas  Santolin  et  al.  [83] 
analysed  the  effects  of  a  real  jet  on  the  cutting  process  of  a 
complete  single-injector  and  rotating  Pelton  turbine,  including 
both  the  penstock  and  casing.  Fig.  11  [83]  shows  details  of  the 
computational  grid  of  the  numerically  analysed  single-jet  Pelton 
turbine  and  of  the  secondary  flows  developing  in  the  real  jet 
outside  the  gate,  generated  by  the  forces  induced  by  the  penstock 
and  the  gate  itself.  The  consequent  deformed  jet  shape  was  shown 
to  modify  the  jet-bucket  interaction  in  comparison  with  an  ideal 
axial-symmetric  jet,  thus  causing  a  low-efficient  energy  exchange 
(hydraulic  efficiency  reduction  of  about  2%). 

Even  though  this  analysis  was  able  to  simulate  an  entire 
working  cycle  of  the  bucket  with  a  good  torque  prediction,  the 
Eulerian  approach,  on  which  it  was  based,  did  not  make  it  possible 
to  analyse  more  in  depth  the  complex  evolution  of  the  free-surface 
flow  pattern  and,  therefore,  appreciate  the  influence  of  the  bucket 
geometry  on  the  turbine  performance.  To  overcome  this  short¬ 
coming,  numerical  simulations  based  on  fully  Lagrangian  meshless 
approaches  were  carried  out  and  attempts  to  simulate  the  flow 
inside  Pelton  buckets  were  also  made  by  adopting  the  Smoothed 


a 


b 


c 


LEADING 

EDGE 


SHROUD 


SHROUD 


Fig.  10.  Comparison  between  numerical  results  obtained  adopting  a  frozen  rotor  interface  model  (a),  a  stage  interface  model  (b)  and  circumferentially  averaged  velocity 
profiles  (c). 
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Fig.  11.  CFD  analysis  of  a  single-jet  Pelton  turbine  [83], 
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Fig.  12.  Results  of  the  hybrid  Eulerian-Lagrangian  approach:  contributions  of  the  different  bucket  discharging  areas  to  the  total  bucket  torque  [87], 


Particle  Hydrodynamics  (SPH)  method  [84,85]  and  the  hybrid 
SPH-ALE  method  [86],  based  on  the  combination  of  the  SPH 
method  and  an  Arbitrary  Lagrange  Euler  (ALE)  one.  However,  even 
though  these  meshless  approaches  have  shown  their  good  cap¬ 
ability  for  predicting  complex  free  surface  flow  fields,  their 
accuracy  still  remains  not  satisfactory  and  further  improvements 
are  required,  as  demonstrated  by  the  lack  of  application  examples. 

An  alternative  and  accurate  method  for  analysing  the  free 
surface  flow  in  a  Pelton  bucket  on  the  basis  of  a  combined 
Eulerian-Lagrangian  approach  was  developed  by  Rossetti  et  al. 
[87],  A  traditional  mesh-based  Eulerian  approach  is  adopted  to 
simulate  the  flow  field  in  the  runner  bucket  in  order  to  exploit  the 
acknowledged  accuracy  and  stability  of  the  mesh-based  methods. 
Subsequently,  to  overcome  the  above  mentioned  analysis  limits  of 
this  method,  the  numerical  results  were  properly  processed  on  the 
basis  of  a  Lagrangian  approach  to  determine  the  particle  trajec¬ 
tories  in  the  runner  bucket  as  well  as  the  variations  of  all  the 
energy  quantities  along  these  trajectories.  This  hybrid  Eulerian- 
Lagrangian  analysis  method  made  it  possible  examine  the  devia¬ 
tion  capability  of  the  different  bucket  areas  and  to  determine  the 
influence  of  the  bucket  geometry  on  turbine  performance  during 
the  jet-bucket  interaction  (Fig.  12)  and  thereby  suggest  possible 
solutions  for  optimizing  bucket  design. 

The  significant  advances  of  the  recent  years  in  the  analysis  of 
the  complex  free  surface  flow  of  Pelton  turbines  have  favored  the 
development  of  innovative  but  undeveloped  design  procedure 
[71,87],  However,  to  define  standardized  and  consolidated  design 


criteria,  as  those  developed  for  other  hydraulic  turbines,  more 
accurate  and  faster  analyses  of  the  complex  free  surface  flow  of 
Pelton  turbines  are  needed  and  a  further  development  of  the  SPH 
method  certainly  represents  one  of  the  most  interesting  challenge. 

4.3.  CFD-based  runner  design  approaches  for  the  rehabilitation 
of  old  small  hydro  power  plants 

In  the  past,  the  design  of  new  turbines  was  usually  based  on 
previous  well  proven,  tried  and  tested  designs  and/or  proper 
model  tests  in  order  to  guarantee  high  efficiency  power  conver¬ 
sion.  In  the  last  years,  the  restrictive  environmental  requirements 
have  significantly  limited  the  development  of  new  plants,  with  a 
consequent  huge  demand  for  the  rehabilitation  and  repowering  of 
old  small  and  micro-hydro  plants  to  satisfy  the  corresponding 
increase  in  the  demand  of  renewable  electricity  production. 
However,  in  the  rehabilitation  and  upgrading  of  existent  old 
plants,  a  new  challenge  had  to  be  faced  as  the  design  of  new 
runners  had  to  match  the  already  existing  stationary  parts. 
Because  of  the  numerous  constrains  imposed  by  the  existing 
components  and  the  uncertainties  connected  to  the  attainment 
of  the  demand  in  power  increase,  the  consolidated,  at  the  time, 
experience-based  design  methods  were  abandoned  and  new  CFD- 
based  design  approaches  started  to  develop.  These  approaches 
swiftly  demonstrated  their  ability  in  analysing  different  and  also 
complex  design  configurations,  making  way  for  significant 
improvements  in  turbine  design.  However,  in  spite  of  the 
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increased  understanding  of  turbine  inner  fluid-dynamics,  the  lack 
of  information  on  how  to  modify  the  runner  blade  for  improving 
turbine  performance  forced  the  designers  towards  time- 
consuming  trial  and  error  approaches,  whose  operating  field  was 
limited  within  the  bounds  of  the  previous  successful  design 
solutions. 

A  significant  improvement  in  the  runner  design  was  achieved 
by  the  development  of  the  3D  inverse  design  methods,  whose 
main  peculiarity  is  to  use  hydrodynamic  parameters,  such  as  the 
blade  loading  [88],  to  determine  the  blade  shape. 

Briefly,  the  blade  loading,  selected  as  the  main  design  para¬ 
meter,  is  fixed  by  imposing  the  distribution  of  the  derivative 
d(rVu)ldm  along  the  mean  meridional  blade  section,  that  could  be 
associated  to  the  pressure  difference  across  the  blade  by  the 
following  expression: 


p+  -p 


2n 

=  pzV" 


d(rVu) 

dm 
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where  p+  and  p~  represent  the  static  pressure  on  the  blade 
pressure  and  suction  side  respectively,  z  is  the  number  of  blades,  p 
is  the  density,  Vm  and  Vu  the  pitch-wise  averaged  meridional  and 
tangential  velocity  respectively. 

Once  Vu  has  been  evaluated  by  the  following  relation: 
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can  be  determined  along  the  considered  blade  section,  where  L/s/U 
is  the  dimensionless  slip  velocity  factor. 

Finally,  for  a  given  radial  distance  from  the  runner/rotor  axis, 
the  angular  coordinate  9  of  the  blade  is  obtained  by 


8  = 


f 


dm 

r  tan  p 
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where  r ^  is  the  radial  distance  of  the  leading  edge. 

Besides  the  blade  loading  distribution  at  two  or  more  spanwise 
sections,  this  method  requires  input  design  parameters:  meridio¬ 
nal  channel  shape,  number  of  blades,  inlet  flow  conditions  in 
terms  of  total  pressure  and  velocity  components,  inlet  and  exit  rVu 
distribution,  stacking  condition  (the  blade  leading  edge  is  not 
aligned  along  the  axial  direction  at  the  rotor  inlet).  These  input 
parameters  perfectly  comply  with  the  design  constraints  imposed 
by  the  existing  stationary  parts  in  the  rehabilitation  and  upgrading 
of  old  plants.  This  peculiarity,  shared  by  all  3D  inverse  design 
methods,  has  significantly  improved  the  effectiveness  of  CFD- 
based  design  methods,  thereby  reducing  the  empirical  aspects  of 
the  former  trial  and  error  methods  and  increasing  the  runner 
design  space. 

The  application  of  these  methods  in  the  runner  design  has 
already  determined  important  design  breakthroughs  such  as  the 
improvements  of  cavitation  performance  and  hydraulic  efficiency 
[88-90]  of  pump-turbine  and  Francis  turbines,  swiftly  demon¬ 
strating  to  be  a  proper  and  revolutionary  design  strategy  in  the 
rehabilitation  of  small  hydroelectric  plants.  However,  further 
improvements  in  runner  design  criteria  are  still  required  in  order 
to  reach  high  efficiency  values  and  to  match  innovative  design 
solutions,  such  as  variable-speed  runners,  with  the  existing 
stationary  parts  of  the  plant. 


5.  Conclusions 

One  of  the  main  goal  of  the  Climate  and  Energy  Policy  of  the 
European  Union  is  to  achieve  a  20%  share  of  renewable  energies  in 
the  overall  EU  energy  consumption  by  2020,  and  the  possibility  of 
increasing  this  share  up  to  100%  by  2050  is  under  consideration  of 
the  European  Renewable  Energy  Council  (EREC). 

In  such  a  context,  hydropower  plays  a  key  role  not  only  as  a 
renewable  and  sustainable  energy  source,  whose  potential  is 
estimated  to  be  still  very  large,  but  also  as  a  large-scale  energy 
storage  technology,  commonly  recognized  as  being  one  of  most 
cost-efficient  among  those  technologies  currently  available.  As  a 
result,  a  renewed  interest  in  pumped  hydro  energy  storage  plants 
and  a  huge  demand  for  the  rehabilitation  and  repowering  of  old 
mini  and  micro-hydro  plants  are  emerging  globally. 

As  regards  pumped  hydro  energy  storage  plants,  the  need  for 
integrating  a  significant  share  of  intermittent  renewable  energy 
sources  in  the  grid  has  turned  the  attention  of  research  towards 
the  development  of  PHES  plants  which  are  as  flexible  as  possible 
in  order  to  provide  frequency  regulation.  To  tackle  this  challenge, 
several  studies  have  been  carried  out  on  innovative  design  solu¬ 
tions  and  advances  in  technology  have  been  obtained  with  the 
development  and  application  of  variable-speed  pump-turbine 
units,  whose  unstable  behaviours  in  both  the  operating  modes 
are  still  under  study.  Moreover,  as  a  consequence  of  the  new  role 
of  PHES  in  the  grid,  optimal  operating  strategies,  based  on  both 
simple  algorithms  and  more  complex  optimization  theories,  have 
been  developed  to  maximize  the  plant  profitability  in  the  deregu¬ 
lated  energy  market. 

The  increasing  need  for  a  better  and  wider  exploitation  of  the 
hydropower  source  has  also  turned  the  attention  of  the  research 
community  towards  the  development  of  more  or  less  complex 
algorithms  for  an  optimal  preliminary  sizing  of  the  plant  capacity, 
on  the  basis  of  economic,  technical  or  techno-economic  para¬ 
meters.  Moreover,  to  increase  the  economic  appeal  of  small-hydro 
power  plants,  and  hence  to  favour  their  further  deployment, 
innovative  management  strategies  have  been  proposed  to  increase 
plant  life  time  and  operating  range,  and  for  maximizing  the 
profitability  when  significant  inflow  and  market  fluctuations 
occur.  A  significant  reduction  in  investment  costs  also  resulted 
from  the  rehabilitation  and  upgrading  of  existent  old  plants.  To 
face  the  numerous  constrains  imposed  by  existing  plants,  innova¬ 
tive  design  approaches,  based  on  computational  fluid  dynamics, 
were  proposed  over  the  last  few  decades,  and  highly  advanced 
design  solutions  were  obtained  in  terms  of  operating  range 
extension  and  hydraulic  efficiency. 

Hydropower  has  all  the  features  needed  for  being  a  sustainable, 
secure  and  competitive  energy  and  for  supporting  a  significant 
penetration  of  other  important  renewable  energy  sources,  such  as 
wind  and  photovoltaic.  However,  even  though  its  future  not  only 
in  Europe  but  all  around  the  world  could  be  bright,  there  is  a  lot  to 
do  to  make  the  technology  more  efficient,  reduce  its  environ¬ 
mental  impact,  modernize  old  plants,  integrate  smaller  plants  into 
the  grid  and  optimize  sizing  and  management.  This  paper  gives  a 
picture  of  the  issues  and  challenges  facing  the  hydro  technology 
sector  and  assesses  the  priorities  for  future  research  and  develop¬ 
ment  with  a  particular  focus  on  pumped-hydro  energy  storage 
plants  and  small  hydro  power  plants. 
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